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Abstract 

X-ray background surveys indicate the likely presence of diffuse warm gas in 
the Local Super Cluster (LSC), in accord with expectations from hydro dynamical 
simulations. We assess several other manifestations of warm LSC gas; these in- 
clude anisotropy in the spatial pattern of cluster Sunyaev-Zeldovich (S-Z) counts, 
its impact on the CMB temperature power spectrum at the lowest multipoles, and 
implications on measurements of the S-Z effect in and around the Virgo cluster. 
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1 Introduction 



As much as 30—40% of all baryons are believed to be in warm gas in large scale 
filamentary structures connecting clusters and superclusters (SCs) of galax- 
ies. Observational evidence for the so-called warm-hot intergalactic medium 
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(WHIM) is still quite rudimentary. First systematic attempts to detect the 
low X-ray surface brightness enhancement that could possibly be expected 
from the WHIM were conducted by Persic et al. (1988, 1990) who searched 
the HEAO-1 A2 database in directions to a sample of SCs. In both these anal- 
yses, as well as a similar analysis of Ginga observations of the Coma-Al367 
SC (Tawara ct al. 1993), only upper limits were obtained on the flux from 
the putative intra-SC (ISC) gas. More recently, evidence for soft X-ray ex- 
cess emission in regions around several clusters has been claimed by Kaastra 
et al. (2003), who carried out a detailed analysis of XMM-Newton line and 
continuum measurements of extended regions around several clusters. 

X-ray emission may be expected, of course, also from the WHIM in the Lo- 
cal SC (LSC); analysis of the HEAO-1 A2 database led Boughn (1999) to 
conclude that this emission was actually detected. Adopting a simple model 
for the morphology of the LSC, and assuming uniform temperature and den- 
sity LSC gas, Boughn deduced the mean electron density to be He = 2.5 • 
10-6 (a/20 Mpcy^^^ {kTe/10 keV)~^'^ cm , when using the specified scalings 
for the values of the electron temperature, kT^, and the length of the semi- 
major axis of the LSC, a. More generally, Kaastra et al. (2003) have reported 
the detection of line and continuum WHIM emission in XMM-Newton ob- 
servations of extended regions around five cluster. While the observational 
results on the properties of LSC g quite uncertain, gaseous filamentary 

structures seem to be a ubiquitous feature of the large scale mass distribution 
described by hydrodynamical simulations (e.g. Jenkins et al. 1998). Of con- 
siderable interest are the statistics and morphologies of these filaments; these 
are currently being quantified (e.g. Colberg et al. 2004). 

The scant observational information on the WHIM and its properties in the 
LSC provides strong motivation for assessing the feasibility of probing it by 
measurements of the Sunyaev-Zeldovich (S-Z) effect (Sunyaev & Zeldovich 
1972). The effect - a spectral change of the Planck spectrum due to Compton 
scattering of the cosmic microwave background (CMB) by electrons in clusters 
of galaxies - is a major cluster and cosmological probe (as reviewed by Rephaeli 
1995, Birkinshaw 1999, and by Carlstrom et al. 2002). In fact, the possibility 
of detectable S-Z signals in directions to SCs has been considered long ago 
(Persic, Rephaeh k Boldt 1988, Rephaeh & Persic 1992, Rephaeli 1993), but 
given the substantial uncertainty in the integrated pressure of ISC, no definite 
predictions could be made. Measurement capabilities have greatly advanced 
since then, so the prospects of detecting weak S-Z signals have correspondingly 
improved. 

Interest in the impact of ISC gas on the CMB has also increased. In particu- 
lar, it has been suggested that gas in the Local Group may be responsible for 
the power suppression of primary CMB temperature anisotropy on large an- 
gular scales. This suppression was deduced from analyses of the COBE/DMR 
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(Hinshaw ct al. 1996) and WMAP databases (Bennett et al. 2003). However, 
the estimated low Comptonization parameter of LG gas makes this suggestion 
quite unhkely (Rasmussen & Pedersen 2001). More recently, Abramo & Sodre 
(2003) proposed that warm LSC gas may be responsible for the deduced CMB 
low multipole power suppression, presumably caused by coincidental align- 
ment of a hot spot in the CMB with the line of sight connecting the Galaxy 
with the Virgo cluster. At the relevant spectral range of the DMR and WMAP 
experiments (~ 20 — 90 GHz), the S-Z diminution lowers the power in the low 
multipoles. Clearly, for this to occur the integrated electron pressure has to 
be sufficiently high. 

WHIM in the LSC may significantly screen our view of the CMB and S-Z sky. 
In addition to the possible impact on the power spectrum of CMB primary 
anisotropy, it will also affect the ability to measure the S-Z induced anisotropy 
and cluster number counts. WHIM in the LSC may appreciably increase S-Z 
cluster number counts. Moreover, the marked ellipsoidal shape of the LSC, 
and the far off-center position of the Galaxy, may result in a substantially 
anisotropic distribution of S-Z counts across the sky. This anisotropy in S-Z 
counts is obviously in addition to that generated by our motion in the CMB 
frame (the CMB kinematic dipole) whose effect on cluster counts was recently 
explored by Chluba et al. (2004). 

This paper is arranged as follows: the LSC gas model, the method for cal- 
culating the directional distribution of S-Z cluster counts, and the relevant 
expressions for calculating the S-Z angular power spectrum, are described in 
section §2. The main results of the calculations are presented in §3. Addi- 
tional aspects of the LSC gas model pertaining to S-Z measurements towards 
the Virgo cluster are discussed in §4. A general discussion follows in §5. 



2 Model and Method of Calculation 

The physical properties of ISC have not yet been well determined; what can 
be done at present is an attempt to incorporate the established observational 
results in a reasonable model in order to provide useful insight to guide upcom- 
ing observations, mostly of the S-Z effect. Relevant LSC properties that seem 
to have been roughly determined are the total mass, baryonic mass fraction, 
and its overall configuration. Knowing the baryon fraction in clusters (e.g. 
Carlstrom et al. 2002), and adopting the scahng to obtain the global baryonic 
fraction (based on hydrodynamical simulations; Evrard 1997), ~ 10 — 14%, 
the total LSC mass, ~ 10^^ Mq, then yields an estimate for the LSC WHIM 
mass, ~ 3 — 5.6 • 10^^ M©, assuming (based on current expectations) that it 
comprises some 30 — 40% of the total mass. The distribution of galaxies and 
clusters in the LSC can be described by a spheroid of semi-axes measuring 



3 



A = 20 Mpc, B = 6.7 Mpc, C = 3.3 Mpc (TuUy 1982); the Galaxy is on the 
major axis of the LSC at a distance of 15 Mpc from the center of the spheroid, 
where the Virgo cluster is located. 

Gas in the LSC is likely to be differently distributed than the galaxies, judging 
by the irregular filamentary structures seen in hydrodynamical simulations. 
However, it is reasonable to expect that most of the WHIM is concentrated 
along the major axis of the LSC mass distribution. Our specific model for gas in 
the LSC is based on the findings of Colberg et al. (2004), who analyzed results 
of N-body simulations (by Kauffmann et al. 1999) in order to characterize 
morphologically the filaments observed in the simulations. They found that 
most filaments have cross sections of 1 — 2 Mpc, with the gas density decreasing 
outward from the major axis of the filament as ~ r^^ beyond some scale radius 
Tg. The LSC gas is assumed to have an ellipsoidal configuration with semi-axes 
(in Mpc) of either 20 x 1 x 1 or 20 x 2 x 2. The highly prolate eUipsoid represents 
the gas 'filamentary' structure. The gas mean density in these structures are 
He = 2.4 ■ 10^^ cm^^, or a factor 4 lower, respectively; its temperature is scaled 
to a value of 1 keV). 

Clearly, the expected anisotropic morphology of ISC gas in the LSC will be 
refiected in anisotropic cluster counts when S-Z surveys are conducted in dif- 
ferent sky directions. The Comptonization parameter towards the Virgo clus- 
ter, a direction which in spherical coordinates corresponds to angles (0, 9) = 
(0,7r/2),is 

y- /aT^ned^ = 3.4.io-« fTinyV 

J meC^ \ \ keV J 



(1) 



whereas the corresponding value calculated for the opposite direction (tt, 7r/2) 
is seven times lower. The non-negligible level of the Comptonization param- 
eter of LSC gas may modify the observed distribution of S-Z clusters across 
the sky; the effective y-parameter measured along a los to a cluster will have 
contributions from both intracluster (IC) and ISC gas, giving rise to an in- 
tensification of the apparent S-Z signal due to the cluster alone. In spherical 
coordinates, the path length through the ISC gas configuration as a function 
of position angle fi = (0, 9) is 
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The Comptonization parameter measured along this direction is then simply 
y{^,9)^aT'^-nsc-r{<P,9). (3) 
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where Tsc and nsc denote the uniform temperature and electron density of 
the LSC gas. 

The flux received from a cluster lying in the direction (0, 6) is a sum of the 
cluster flux and the S-Z signal generated in the LSC; assuming a /3-King 
profile for the IC gas density with /3 = 2/3 and an isothermal temperature 
distribution, the total flux at a frequency v is 



AF((/), 6) = AFc{(P, 6) + AF5c(0, 9) = 



2 {kT^f / k(jT ^ 
{hcf ' \ ~^\ 



g{x) 
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where Tcmb = 2.726°-?^, Tq, Uq, Tsc ^-nd nsc denote the temperature and den- 
sity of the cluster and LSC gas, respectively, p = Rmr/f^c, <^b is the beam width, 
and g{x) is the non-relativistic spectral form of the thermal S-Z effect, valid 
here owing to the relatively low temperatures of both the ISC and Virgo IC 
gas. At higher temperatures a relativistically correct expression (e.g. Shimon 
& Rephaeli 2004) must be taken into account. 

The calculations were carried out in a standard ACDM model, specified by 
the parameters — 0.3, 0,a — 0.7, h — 0.7, n — l,as — 1. We use the cluster 
mass function of Sheth & Tormen (2001) mass, and the temperature-mass 
relation 



where E"^ = + z)^ + ^a- A non-evolving gas fraction /q = 0.1 is as- 

sumed. And since our treatment here is most relevant for the Planck all sky 
survey, we have selected the spectral band to be the HFI uo—353 GHz, with 
Au — 116 GHz; correspondingly, the beam profile is taken to be Gaussian 
with FWHM of 7.1'. 



Number counts are calculated based on the expression 

,dr 

d; 



N{> AF^) = j r^^dz j N{M,z)dM, (6) 
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where the lower hmit is assigned such that the flux received from a cluster with 
mass M situated at redshift z exceeds the flux limit of the survey. In order 
to estimate the number of clusters that may be observed above a flux limit 
within the framework of the LSC model, the LSC S-Z flux is mapped using the 
second term in equation (4); the results are illustrated in the upper left- (model 
1) and right-hand (model 2) panels of figure 1. The number of clusters with 
flux exceeding the limit AF^, are then calculated using equation (6). Calculated 
number counts are then multiplied by ~ 6-10"*^ to yield the number of clusters 
expected within square patches of sky measuring 5° x 5°. These are depicted 
in the lower panel of figure 1. For the detection of a cluster with a given fiux to 
be possible, the combined S-Z fiux from the cluster and the intervening LSC 
gas must exceed that of the measurement sensitivity. For example, if the latter 
is 30 mJy, then in a sky region where the S-Z fiux due to the LSC is 20 mJy, 
the minimum detection level for a cluster is 10 mJy. Having determined the 
minimum fiux, the expected cluster counts within a 25 square degree patch of 
sky in the requested direction can be read off the graph in the lower panel of 
figure (1). 

The angular power spectrum of the S-Z effect due to LSC gas can be calculated 
using the usual spherical harmonic expansion of the temperature anisotropy. 
Here, however, the flat sky approximation is invalid due to the large angular 
scales involved. The basic expression is 



which, after projecting the three dimensional function AT/T(r, 0, 6) onto the 
two dimensional celestial sphere, and some algebraic manipulation, assumes 
the form 



where r(0, 9) is given by equation (2). The angular power spectrum may then 
be easily calculated as 
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Fig. 1. Upper left- and right-hand panels: color-coded flux levels of the S-Z signal 
at 353 GHz (in units of mJy) from LSC gas at a density He = 2.4 • 10~^ cm~^, and 
a density which is a factor 4 lower, respectively. Lower panel: number of clusters 
within 5° x 5° patches of sky, above flux limit AF,^. 

3 Results 

For flux limits of 30 mJy and 60 mJy the corresponding all-sky cluster counts 
sum up to ~ 40 000 and ~ 4 300, respectively. Assuming that the cluster 
population is homogeneously distributed across the sky, this amounts to ~ 24 
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and ~ 3 clusters in each 25 square degrees (rectangular sky region). It is 
important to emphasize that the number of potentially detected clusters in 
future S-Z surveys predicted by numerical studies is quite sensitive to the 
choice of parameters characterizing the background cosmology and cluster 
properties, and consequently, these counts merely reflect the specific modeling 
detailed in the last section. Now suppose that observations are conducted in 
the directions ~ (0, 9) = (0°, 90°) (i.e. towards the longest possible los within 
the LSC halo), (5°, 90°), (180°, 90°) (designating the direction opposite to the 
Virgo cluster), and, e.g., (180°, 20°). The S-Z fluxes induced by the LSC gas 
in these directions are (see figure 1) ~ 21, 7, 3, O.BmJy, respectively, implying 
that for a flux detection limit of 30 mJy clusters will be detected if their S-Z 
signals exceed the complementary values (9, 23, 27, 29.5 mJy, respectively) in 
these directions. The corresponding limiting values for a mean gas density of 
rie = 6 • 10-^ cm-^ are ~ 25, 27, 29.3, 29.7 mJy, respectively. 



The number of clusters that could be detected in these directions, i.e., those 
which meet the detectability criterion, can be easily read ofi^ the diagram in 
figure 1. For example, an experiment with a detection limit of 30mJy yields 
~ 570, 73, 50,40 clusters per 25 square degrees. The corresponding counts for 
the larger ellipsoid (with gas density lower by a factor of 4) are ~ 60, 50, 40, 40. 
A limiting fiux of 60 mJy yields counts of ~ 20,8,7,6 in the first case, and 
~ 6 — 8 for all directions in the second case. While these estimates of the clus- 
ter number counts are substantially uncertain, and the fact that there is also 
considerable variation along difi'erent los across 25 square degrees (particularly 
towards both ends of the spheroid major axis), the results nonetheless demon- 
strate the possible impact of LSC gas, especially in a survey with 30 mJy fiux 
limit, and in direction to the Virgo cluster. The difference between counts 
along los in the Virgo region and those in other directions markedly exceed 
the statistical uncertainty (of a purely Poissonian distribution of clusters) . 



The angular power spectrum of the S-Z effect induced by LSC gas is illustrated 
in figure 2. The maximum power level is significantly lower than the corre- 
sponding values found by Abramo & Sodre. This is due to our lower values 
of the density and temperature, each lower by a factor of ~ 2 than their val- 
ues (5 ■ 10^^ cm^"^, and 2 keV), and the very different LSC gas configuration 
assumed by Abramo & Sodre, which is reflected in a different distribution 
of power among the explored multipole range. In particular, more power is 
seen at higher multipoles (with regard to Abramo & Sodre's results), which 
is expected owing to its highly prolate morphology. Consequently, given the 
substantial uncertainty in the properties of ISC gas, one cannot convincingly 
argue that the observed suppression of the primary CMB temperature power 
spectrum is due to Compton scattering in the LSC. 
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Fig. 2. The angular power spectrum of the S-Z effect induced by the gaseous LSC 
halo. 

4 Implications for the Virgo Cluster 



The central location of Virgo in the LSC imphes that measurement of the S-Z 
effect in this cluster will likely include a substantial contribution from LSC 
gas. The Comptonization parameter measured along the direction (0, 9) is 



(10) 



where r(0, 9) is given in equation (2) and 
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p2 - (9/9,)^ ■ H^F, 



13 -{p'-{9/9^r) 
2' '2' l + {9/9c)^ 



is the integrated line of sight through the Virgo cluster. Here 9c is the core 
radius, 6 = 3/2/?7 where 7 is a polytropic index, p is the virial to core radius 
ratio, and H2F1 is the hypcrgcometric function. The main gas configuration 
in Virgo is centered on the giant elliptical M87, to which we refer here. Ex- 
cluding a relatively small region in the center of M87, the gas is taken to be 
isothermal with a (King) P density profile, having (Matsumoto et al. 2000) 
no = 0.019 cm-^ Tq = 2.28 keV, = 0.014 Mpc, p = 125, and (3 = 0.4. The 
magnitude of y along the central Virgo los increases by ~ 20% due LSC gas, 
its gradient across the cluster is appreciably shallower than would be expected 
from the intrinsic density profile, and the cluster seems slightly larger due to 
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the enhancement of the cluster Comptonization parameter. 



The combined motion of the earth towards Virgo and the motion of the LSC 
in the CMB frame accounts for the dipole term; the deduced velocity is ~ 
600 kms^^ in the direction (£ — 264°, b — 48°) in galactic coordinates (Fixen 
et al. 1996). The temperature change generated by the dipole in the direction 
of Virgo is AT/T ~ 1.8 ■ 10^'^, whereas the combined S-Z temperature change 
due to the LSC halo and Virgo cluster peaks at AT/T ~ 5 • 10~^. Clearly, the 
dipole may interfere with S-Z measurements of the Virgo cluster and should 
therefore be taken into account. This would generally be the case for large 
angular diameter clusters, across which the dipole gradient can be significant. 
The temperature change due to the dipole term is simply 



whereas the temperature change due to the thermal S-Z effect in the R-J 
region is AT/T = s{x) ■ y, where s{x) = x/ tanh(a;/2) — 4. In order to be able 
to compare the contributions of both effects to the temperature change, the 
galactic coordinates specified above are converted into equatorial coordinates, 
and transformed such that the Virgo cluster lies in the direction (0 = 0°,^ = 
90°). In this coordinate system the velocity vector points at ~ (0 = 340°, 9 = 
110°). Figure 3 illustrates the relative temperature change (AT/T) due to 
the LSC-I- Virgo S-Z signal, and the latter combined with the CMB dipole 
{AT/T)dipoie+i^T/T)sc+Virgo withiu a 10° x 10° patch of sky centered around 
the Virgo cluster. The relative temperature change due to the dipole does not 
change considerably across the patch due to its small angular size and the fact 
that the angle formed between the dipole and the vector pointing at Virgo 
lies within the range 21° < 7 < 35°, at a significant angular distance from 
the direction at which the cosine term of the dipole has the steepest slope, 
7 = 90°. The angle between the two vectors can be calculated using 

cos 7 = cos^i COS ^2 + sin^i sin6'2 cos (0i — 02)- (13) 

The combined relative temperature change is still dominated by the dipole 
term, although its profile is now deformed due to the S-Z induced temperature 
change. In particular a 'warm spot' may be seen in the direction of Virgo 
(which translates into an enhancement of ~ 3% of the temperature change due 
to the dipole term in this direction), by virtue of the fact that at a frequency 
of 353 GHz the cluster behaves as a source for CMB photons. Consequently, it 
is obvious that the dipole term must be taken into account in measurements 
of S-Z signals in nearby clusters. 

The impact of LSC gas on measurements of the S-Z towards the central Virgo 
region is appreciable only when observations are made in a survey mode - such 




(12) 
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as planned with PLANCK - with relatively short scan time over the region, 
and consequently insufficient sensitivity to remove such an LSC component. 
Clearly, the additional signal due to LSC gas can largely be subtracted out 
when pointed observations are made with a suitably selected beam-through 
pattern. 



xlO'* x10' 




Fig. 3. The distribution of the relative temperature change around the Virgo cluster 
due to IC and ISC gas (left-hand panel), and combined with the temperature change 
induced by the CMB dipole term (right-hand panel). 



5 Discussion 

Although there is yet no unequivocal observational evidence for the WHIM, 
theory and hydrodynamic cosmological simulations indicate that it is highly 
prevalent and constitutes a substantial fraction of the baryonic mass fraction. 
Currently available observations of WHIM in several cluster regions imply 
appreciable levels of the Comptonization parameter, ~ 10~^, comparable to 
values in most clusters. Even though we used a simplified model for the WHIM 
in the LSC in order to explore some of its consequences, we consider the re- 
sults presented here to be qualitatively valid. In fact, our approximation of 
the WHIM filamentary structure by an ellipsoidal configuration (with a large 
volume) quite likely underestimates its S-Z impact, given the basic premise 
that it constitutes an appreciable fraction of the baryon mass of the LSC. As 
indicated by morphological analyses of such filaments seen in hydrodynamical 
simulations, their density grows towards the central axis, such that a clear 
preferential direction marked by the axis can generate even higher S-Z signals 
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along this direction. Moreover, the adopted symmetric morphology underes- 
timates the density along the main axis in direction to Virgo as there is no 
cluster in the opposite side of the LG- Virgo direction. 

Several potential observational consequences of the presence of such gaseous 
component have been discussed in this paper: the angular power spectrum 
of the induced S-Z signal has been calculated and shown to be rather low, 
particularly with regard to what has been speculated in the literature as a 
possible explanation for the suppression of primary CMB power on the lowest 
multipoles. While we cannot rule out a suppression at the claimed level, we do 
not find this to be very likely. The effect of LSC gas on directional S-Z cluster 
counts is proportional to the los crossing the volume of gas in the LSC. On 
the other hand, a definite prediction can be made that the distribution of S-Z 
cluster counts across the sky will be determined to be anisotropic close to the 
level estimated here. The impact is relatively significant owing to the highly 
non-linear form of the mass function; a relatively low S-Z fiux induced by 
intervening LSC gas may give rise to a considerable increase in the number of 
potentially detected clusters. This is in particular true in light of the decidedly 
higher abundance of low-flux S-Z clusters predicted by the mass function. 
Needless to say, this effect is most noticeable along the main axis towards 
Virgo. S-Z measurements of the Virgo cluster may also be affected by the 
presence of the WHIM, as well as by the dipole distribution of the primary 
CMB temperature anisotropy. This is even more relevant to measurements of 
the kinematic effect, since the induced temperature change does not depend 
on the temperature of the LSC gas, which tends to be lower than typical 
IC temperatures. These effects will have to be taken into account in realistic 
analyses of the S-Z signal from the Virgo cluster. 

Finally, although observational evidence for the WHIM is still scarce, it is 
interesting to use the currently available data to assess its potential of gener- 
ating a significant S-Z signal towards other SCs. Kaastra et al. (2003) list five 
clusters (A1795, Sersic 159-03, MKW3s, A2050, A2199) reported to include a 
WHIM component, modeled to he within a spherical halo, be isothermal and 
of homogeneous density. For each cluster they list the WHIM electron density, 
temperature, and size, such that the Comptonization parameter along the ra- 
dius of the spherical halo can be easily calculated. They all turn out to lie 
within the range y ~ (0.4 — 0.8) • 10~^, i.e. near the level taken here for the 
LSC. Thus, the S-Z consequences of ISC gas we considered here are hkely to 
be relevant in other nearby SCs, but obviously on smaller angular scales. 
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